A=[1 -1 2 -1;-1 1 3 -2;2 3 1 0;-1 -2 0 1];
[V,D]=eig(A);
disp('特征值为：');
disp((diag(D))');
disp('对应的特征向量为：');
disp(V);